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Abstract 

The loop level lepton flavor violating signals 77 — ► ££' (£ = e, fj,, r, £ ^ £') are studied in a scenario 
of low-energy, R-parity conserving, super symmetric see-saw mechanism within the context of a high 
energy photon collider. Lepton flavor violation is due to off diagonal elements in the left s-lepton 
mass matrix induced by renormalization group equations. The average slepton masses fh and 
the off diagonal matrix elements Am are treated as model independent free phenomenological 
parameters in order to discover regions in the parameter space where the signal cross section may 
be observable. At the energies of the 77 option of the future high-energy linear collider the signal 
has a potentially large standard model background, and therefore particular attention is paid to 
the study of kinematical cuts in order to reduce the latter at an acceptable level. We find, for the 
(er) channel, non-negligible fractions of the parameter space (5ll = Am 2 /fh 2 > 1CP 1 ) where the 
statistical significance (SS) is SS > 3. 

PACS numbers: ll.30.Fs, 11.30.Pb, 12.gO.Jv, 14.80.Ly 
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I. INTRODUCTION 



The high-energy linear lepton collider (LC) is presently considered as a necessary next 
step in the field of high-energy physics. If new physics will show up at the CERN Large 
Hadron Collider (LHC), a LC with a much cleaner environment would allow unambiguous 
precision measurements. However the LC project has the potential to address, on its own, 
questions of physics beyond the standard model, since e~e~ and 77 options are also planned 
beside the basic e + e~ mode. If these options are carried on, they will provide us for the 
first time with the high physics potential of very high-energy e~e~ and 77 collisions. See 
for example [l] for a full discussion of the physics potential of the TESLA photon collider 
(PC). 

A topic which has recently received considerable attention is that of neutrino mass and 
lepton number (flavor) violation, LNV (LFV). Non- vanishing neutrino masses induce LFV 
processes such as I — > £'7. If neutrinos have masses in the eV or sub-eV range, the neutrino 
generated branching ratio to the latter process is of order (9(10 -40 ) and therefore unobserv- 
ably small. For such processes to be experimentally accessible, new physics has to come into 
play. Experimental searches of radiative lepton decays put strong bounds on models of LFV: 
Br(n -> erj) < 1.2 x 10" 11 Q, Br(r -> erf) < 3.9 x 10" 7 Br(r -> yrj) < 3.1 x 10~ 7 Q. 
Supersymmetric (SUSY) extensions of the SM in the soft SUSY breaking potential V so f t con- 
tains, in general, non diagonal entries in generation space and therefore additional potential 
sources for LFV. Even in minimal supergravity scenarios characterized by universal soft 
mass term for scalar slepton and squark fields, renormalization induces potentially sizable 
weak scale flavor mixing p| in V so ft- 

In this paper we study the lepton flavor violating reaction 

77 -> t£ (1) 

with i 7^ £' and £, £' = e, /1, r, which arises at one loop order in the just mentioned SUSY 
scenario, thus extending to the 77 option an analysis done by some of the authors in Ref. []| 
for the e + e~ and e~e~ mode of the next linear collider. The OPAL collaboration searched for 
this type of LFV reactions up to the highest center-of-mass (CM) energy reached by LEPII, 
\fs = 209 GeV |f|. One e + e~ — > efi event was found at yfs = 189 GeV matching all tagging 
conditions, but it was interpreted as due to initial state radiation. These processes have 
the advantage of providing a clean final state which is easy to identify experimentally (two 



back-to-back different ffavor leptons), though one has to pay the price of dealing with cross 
sections of order 0(a A ). In Ref. [?J we found that the e~e~ option stands better perspectives 
for the possible detection of a LFV signal as opposed to the e + e~ mode, both because of 
larger cross sections and smaller background. In general the 77 mode offers larger cross 
section as compared to the other modes, but at the same time has the drawback of larger 
background and one must take into account the non-monochromaticity of the beams. 

The plan of the paper is the following: in section |H] we discuss the SUSY scenario of LFV 
in the charged slepton sector (details of the helicity amplitudes of the diagrams contributing 
to the signal reactions are given in the appendix); in section ITTT1 we review briefly the photon 
spectra used in the numerical computations of the signal ; in section HVl we discuss the main 
features of the signal; in section E] we discuss the main standard model (SM) backgrounds, 
and finally in section IVO we present the concluding remarks. 



II. SUSY SCENARIO FOR LEPTON FLAVOR VIOLATION 

In the SUSY extension (with mSUGRA boundary conditions) of the seesaw mechanism 
for the explanation of neutrino masses [8|, the superpotential contains three SU(2)l singlet 
neutrino superfields Ni with the following couplings jfl lil llfl|: 

W = (Y v )ijE a pH£NiLj + ^(M^NiNi. (2) 

Here H 2 is a Higgs doublet superfield, Lj are the SU{2) L doublet lepton superfields, Y u is a 
Yukawa coupling matrix and Mr is the SU(2)l singlet neutrino mass matrix. At low energy 
the renormalization group equations (RGE) produce within the MSSM diagonal slepton 
mass matrices. With the additional Yukawa couplings in Eq. (J2J and the new mass scale 
(Mr) the RGE evolution of the soft SUSY breaking parameters is modified : assuming a 
heavy right handed singlet neutrino mass scale, Mr, the RGE from the GUT scale down 
to Mr induce off-diagonal matrix elements in (rrv~)ij. In the one loop approximation the 
off-diagonal elements are jjj: 



K)y * -^(3 + al)ml(YXh ^ (^gf ) • (3) 



ao is a dimensionless parameter appearing in the matrix of trilinear mass terms A? = Yf aomo 
contained in V so ft- The rate of LFV transitions like ti — > £j, i j, £ — e, fi,T induced by 
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the lepton-slepton-gaugino vertex is determined by the diagonalization matrix Unj. These 
matrix elements can be potentially large because they are not directly related to the mass of 
the light neutrinos, but only through the seesaw relation m v ~ m 2 D /Mf> = v 2 Y 2 /Mr. The 
same effect on the mass matrix of SU(2)l singlet charged sleptons (jriQij is instead much 
smaller : indeed, in the same leading- log approximation of Eq. (j3J, the corresponding RGE 
do not contain terms proportional to YJY V , since the right-handed leptons fields only have 
the Yukawa coupling Yg, which completely determines the Dirac mass of the charged leptons 
and these are known to be small numbers. Thus the off-diagonal elements of (m^)y can 
be taken to a vanishing to a very good degree of accuracy. The mixing matrix arising in 
the diagonalization of induce LFV couplings in the lepton-slepton-gaugino vertices 

i^.UiijiLjX- The magnitude of LFV effects will in turn depend on the fundamental theory in 
which this mechanism is embedded (for example SU(5) or 50(10) SUSY GU T [lol 111 I12I]) 
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and on the particular choice of texture for the neutrinos mass matrix 

In this paper we adopt a more general and phenomenological approach, as it was done 
in jzj], without referring to a particular GUT model or neutrino mass texture. We consider 
a two generation model for the mass matrix of left-sleptons (and sneutrinos): 

/ m 2 Am 2 \ 




' 2 

m 



with eigenvalues: m 2 ± = m 2 ± Am 2 and maximal mixing. Under these assumptions, the 
lepton flavor violating propagator (in momentum space) for a scalar line is 



1 1 \ Am 



,2 



(^j)o 2\p 2 — m 2 + p 2 — m 2 ) \p 2 — m 2 + ){p 2 — m 2 __) ' ^ 



while the propagator for a lepton flavor conserving (LFC) scalar line is: 



(6) 



2 \p 2 — m\ p 2 — m 2 _ 
The quantity 

S LL = Am 2 /fh 2 (7) 

is the dimension-less parameter that controls the magnitude of the LFV effect. This approach 
allows us to study the signal in a quite model-independent way by means of scans in the 
parameter space - the (in, 5ll) plane - which is already constrained by the experimental 
bounds on radiative lepton decay processes. 



For the calculations presented in this work it is a good approximation to assume that the 
two lightest neutralinos are pure Bino and pure Wino with masses Mi and M 2 respectively, 
while charginos are pure charged Winos with mass M2 , Mi and M2 being the gaugino masses 
in the soft breaking potential. The Higgsino contribution to neutralino and charginos has 
suppressed amplitude, since the coupling is proportional to the lepton masses. For the 
same reason left-right mixing in the slepton matrix is neglected. The relevant parts of the 
interaction lagrangian are, adopting the notation of 



£ = -g£P R WD + -^=£P R W 3 £ L + -^=t w £P R B£ L + e£ L A^L*i d H L - eA^W-fW ( 8 ) 



The contributing one loop diagrams are displayed in Figure ^ We have grouped them 
according to their topology : (a) penguin type ; (b) self-energy ; (d) box diagrams. It is 
of course understood that each diagram is accompanied by an exchange diagram where the 
final state leptons (or initial state photons) are exchanged. 

The possibility of having, at the next LC, high-energy polarized photon beams suggests 
(see discussion in section IIII|) to calculate the amplitudes of the diagrams in Fig. ^ within 
the helicity formalism. Denoting by T, X, Y, Z the space-time unit four-vectors, the four- 
momenta of the particles in the center-of-mass frame (CMF) are expressed as: 



Pl= 2 (^ + ^) P3 = } y(t + X smO* + Z cos6* 
p 2 = xl(f-z) p 4 = ^ (f-Xsine* -Zcos9* 



(9) 



where s = (pi + P2) 2 arid e* are, respectivley, the CMF energy and scattering angle, while 
the polarization four-vectors of the photons are: 

e > = _-L(Ax + ;y) 4 = +-j= (x'Jt - if) , (10) 

where A and A' (= ±1) denote the photon helicities. Assuming massless external fermions 
and, given the chiral nature of the coupling in the lagrangian, the helicity of the fermions in 
the final state are fixed to only one configuration, thus there are only four helicity amplitudes 
corresponding to the possible combinations of the photon helicities. With obvious notation, 
we indicate with M.( x ' x '^ the helicity amplitudes. The loop integrals are decomposed in form 
factors according to the notations of the software package LoopTools [l^ which is used 
in the code for numerical computations. The final analytical formulas of the amplitudes 
j^(x,x') are f unc tion of s, 9*, A, A' and the SUSY parameters. They also contain form 



factors originating from the loop integrals which are defined according to the LoopTools 
conventions. We report the explicit expressions of the helicity amplitudes in the Appendix 1X1 





FIG. 1: Diagrams for 77 collisions : (a) penguin diagrams ; (b) Self-energy diagrams ; (c) box 
diagrams. The full black dot in a scalar line denotes the insertion of the lepton flavor violating 
propagator (Eq. 0). In the diagrams of part (c) this insertion is to be done successively in each 
scalar line. In addition all graphs are accompanied by an exchange graph where the final leptons 
are interchanged. 
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III. DISCUSSION OF PHOTON BEAMS AND PC LUMINOSITY 



High-energy photons beams |18| will be obtained from Compton back-scattered (CB) 
low-energy laser photons with energy uq off high-energy electron beams with energy Eq. 
These high-energy photon beams will not be monochromatic but will present instead an 
energy spectrum, mainly determined by the Compton cross section, up to a maximum 
energy y m E , where y m = x/(x + 1) with x = AEquiq/w^. 

Full simulations of the experimental apparatus, see for example the simulation of Telnov 
for TESLA fl, show that the real luminosity spectrum cannot be described by simple 
analytical formulas because of energy-angle correlation in Compton scattering, collisions 
effects and details of the collision region. Besides the high-energy peak also a 5-8 times 
higher low-energy peak is present, which is originated by photons after multiple Compton 
scattering and beamstrahlung that cannot be described by analytical formulas. 

The high-energy peak is instead found to be almost independent of the technological de- 
tails and well reproduced by the product of two Compton spectra. The normalized Compton 
energy spectrum is: 



Fix y) = ^ d -^ = ^ 



' y + {2r-lf -\ e P i xr{2r-l){2-y) 



where N c is the normalization constant 1 , y — E 1 /Eq is the fraction of the initial electron 
energy acquired by the CB photons, r = y/x(l — y), and A e and Pi are the electrons and laser 
photons polarizations (|A e , Pj| < 1), respectively. Thus the theoretical differential spectrum 



for luminosity is: 



dL CB 

- F c {x, yi )F c {x,y 2 ). (12) 



dy 1 dy 2 

It is useful to rewrite Eq. (fT2"j) in terms of the invariant variables z = ^/yvy 2 = W 11 /2E 



\J s^-y/sgg and the pseudo-rapidity r] = In ^/yi/y 2 , and define a differential spectrum as a 
function of z : 

dL CB rlny m /z 

—^ = 2z F c [x,ze^)F c {x,ze-^)di]. (13) 

">Z J-\ny m /z 

This is the function we have plotted in Fig. 121 for some values of Eq and for correlated value 



1 Nr. 



I" I -Jr)Mz + l) + i + i- +A e P i [(l + i)ln(* + l)- § + ^ 



l+x 2(> + l) 2 
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FIG. 2: The ideal CB luminosity spectrum plotted for different CM energies of the ee collider. 

of x calculated using the TESLA parameter ujq = 1.17 eV. It gives a peak of luminosity 
near the maximum value of z, z m = W limax /2E = y m , as shown in Fig. (J2J) and a broad 
spectrum at lower values. This means that most of the collisions involve two high-energy 
photons from the high-energy peaks of the two Compton spectra. We note that the peak in 
the luminosity spectrum is obtained when the product Pi X e is negative for both Compton 
spectra. 

In this high-energy range, the colliding photons have practically the same energy which 
is close to its maximum value. Obviously, this configuration is the most favourable to 
distinguish two-particle final states among multi-particle production. It is important to 
notice that the experimental design for a future photon-photon collider is planned so as to 
have full control of the luminosity and optimize it in this high-energy range in view of Higgs 
physics studies |19]. In the low-energy range, collisions between photons that may have 
very different energies take place, leading to copious boosted events. Then, the separation 
between signal and background becomes more challenging. Moreover, this low-energy part 
is more dependent of the experimental apparatus. For these reasons, we have restricted our 
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study to the high-energy part of the luminosity spectrum. Another reason to restrict the 
peak is that the total luminosity of the photon collider is defined by the condition 

L 77 = / dz^TL. (14) 

To evaluate the expected total number of events and event rates we take as benchmark 
the TESLA parameters in Ref. |l|. At \f^ee — 2E = 200, 500 GeV the geometrical lu- 
minosities are expected to be L = (4.8, 12, 19.1) xlO 34 cm -2 s -1 , while the corresponding 
photon-photon luminosities at the peak are : L 77 (z > 0.8z m ) =(0.44, 1.15) xlO 34 cm~ 2 s _1 , 
equivalent to (1.3, 3.4) x 10 2 fb _1 yr _1 . To use these simulated realistic numbers with the 
ideal spectrum in Eq. (fTT?|) a suitable normalization is necessary : 

fZmax Any m /z 

L 77 = Cnorm / 2z dz / F c (x, ze +v )F c (x, ze~ v )dr]. (15) 

Jo.Szmax J — \ay m /z 

pZmax dL 

N events = L 77 / dz- v I1 o-(W 77 ) (16) 

J0.8z max L n az 



Substituting Eqs. ()15pi6j) into the integral we eliminate the dependence from C norm , which 
depends on the total integrated luminosity, redefining the differential spectrum as 

J T norm J j CB 

dz norm dz ' 1 } 

where L norm is given by: 

Thus we define both for signal and background the effective cross section as : 

^.effective _ / dz _T^ (y{ y V ^ ( 19 ) 

and the total number of events is thus given by N events = £ 77 x <j e ff ecttve . In view of studying 

helicity correlations, we discuss the polarizations properties of the back-scattered photons. 

The degree of circular polarization is given by: 

= -P t (2r - 1) [(1 - y)- 1 + l-y]+ 2\ e xr [1 + (1 - y)(2r - l) 2 ] 
{ l} N c F c (x,y) 1 ' 

Assuming complete polarization for laser photons (Pi = ±1) and the planned maximum 

available for electrons A e = ±0.85, this function is plotted in Figure El for various values of 

x. As can be seen, in the high-energy peak where y is near y m , colliding photons have a 

high degree of circular polarization with P 7 = —Pi. 
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FIG. 3: The ideal CB helicity spectrum plotted for different center of mass energies of the ee 
collider. 

IV. DISCUSSION OF THE SIGNAL 



We discuss first the signal for the ideal case with monochromatic photons in pure helicity 
state. The differential polarized cross sections with respect to the scattering angle in the 
photon-photon CMF are given by 

da xy 1 



\M XX ' (s,t,u)\ 2 



(21) 



dcos9* 327TS 1 

We plot them in Figure 0] as functions of the CMF scattering angle with masses set to the 
values specified in the caption of the figure and for x /sZZ =128 GeV which corresponds to 



77 

the maximum energy that is reachable with a LC with y/s^ = 200 GeV. It is seen that 
the amplitudes with opposite helicity photons M.^ + ~^ and (J z = ±2) dominate 

the signal, while those with same helicity photons (J z = 0) give negligible cross-sections. 
Moreover the former are peaked in the forward and backward directions while the second are 
suppressed in these regions. The total cross sections are plotted in Figure El varying the CM 
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COS0 



FIG. 4: Differential cross section given by the four helicity amplitudes for monochromatic photons 



at 



'77 



128 GeV. The values of the masses are Mi = 200, M 2 = 100, (m f ) = 150 GeV and 



Am 2 = 6000 GeV 2 



energy. The J z = ±2 cross sections decrease with increasing energy for they are dominated 
by the diagrams with the exchange of a light lepton in the t and u channels [(a) and (b) of 
Fig. P. 

The realistic effective differential cross sections as function of the scattering angle in the 
laboratory (e~e~ CMF) are simply obtained by a boost and by convoluting the "monochro- 
matic" differential cross-sections in Eq. (f2~T]) with the luminosity spectrum discussed in the 
preceding section. The fact that photons are not in pure helicity state is here taken into ac- 
count using density matrices for initial photons expressed in terms of Stokes parameters (^(J . 
The complete formula is: 

d° xx ' [ Zma * ^™rfcosr(l-(A(x,P, 1 ,A ei )))(l-(A'(x,P, 2 ,A e2 ))) da^' 

dz JZ TZZT a o n ( 22 ) 



dcosO } v dz dcos9 2 2 dcosO* 

Here the functions (A(x, P^, A ei )), (A'(x, Pi 2 , A e2 )) play the role of the Stoke parameter r] 2 , 
while rji and r/3 [20] give no contribution for we are assuming laser photons with full circular 
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FIG. 5: Total signal cross sections for monochromatic photons as a function of the energy with the 
parameters as specified in the caption of Figure 0] . 

polarization. The total cross sections are obtained finally by integrating over the laboratory 
scattering angle and introducing a kinematical cut : 



In Figure El we study the effect of the inclusion of spectra in the calculation. The upper 
panel presents the same monochromatic (+ -) differential cross section of Figure 0] com- 
pared with the one calculated with the complete formula of Eq. E21 while the bottom panel 
shows the total cross section as a function of the parameter 5ll (c.f. Eq. |7j) : it is clear 
that the complete formula in Eq. E21 gives results almost identical (within a few %) to the 
monochromatic calculation with photon energies fixed to their maximum values. This is 
a consequence of the choice of restricting the calculation to the luminosity peak near z m 
where the photons have energies near E™ ax , are in an almost denned helicity state and the 
boost to the lab system has (3 ~ 0. Thus, in the following we consider the cases of a photon 




(23) 
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FIG. 6: Effect of energy and helicity spectra on angular and total cross sections for the helicity 
configuration (+ -). The values of masses are the same as in Figures HJ 03 



collider with 2Eq = 200 and 500 GeV, with monochromatic photons in pure helicity states 
with E 1 = E™ ax (y/s^ = 128 and 410 GeV respectively) and use the realistic simulated 
luminosities of TESLA to estimate event rates. 

In Figures [7| and |H1 we plot the cross section given by the dominant amplitude (+ -) as 
function of the insertion 5ll for some values of gaugino and average slepton masses. These 
values are in the range of the SPS1 benchmark point mSUGRA scenario 21\ that give a 
particle spectrum with the lightest charginos, neutralinos and sleptons in the 100 — 200 
GeV region. This light spectrum is also favoured by global fits to the standard model 
parameters j^J. Even if the differential cross section is peaked along the collision axis, a 
necessary angular cut | cos#| < 0.9 is applied because the background is also large in this 
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(M 2 , Mj)=(180, 90) GeV (M y Mj)=(200, 100) GeV 




(M 2 , M 1 )=(300, 150) GeV (M y M 1 )=(400, 200) GeV 




FIG. 7: Total cross section for the amplitude (+ -) as a function of the dimensionless parameter 
5ll and ^/s 77 = 128 GeV. The values of the others parameters are given in the legends. 



region, as it is discussed in Section IVl Given luminosities of order (9(100) fb _1 yr _1 , cross 
sections greater than 10~ 2 fb are needed. In the case of a 2Eq = 200 GeV PC, this happens 
for Sll > 5 x 10~ 2 while in the 2E = 500 GeV case, 5ll > 5 x lO^ 1 is needed, which means 
quite large non-diagonal matrix elements. 

To see if these large mass splittings are allowed by current experimental constraints 
we have to take into account the bounds imposed on the model by the non observation 
of radiative decays. In Figures |Hl EH we show scatter plots where the average lepton 
masses and the relative mass splitting 5ll = Am 2 /rh 2 are varied freely, for fixed values 
of gaugino masses. All the parameter space [the (fh, plane] is covered by the clear 
triangle- up shaped points (turquoise in color) that satisfy the bounds T(r — > 7//, (e)) < 
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(M 2 , M 1 )=(180, 90) GeV (M y M 1 )=(200, 100) GeV 
!| I I 1 | I I 1 




(M 2 , M 1 )=(300, 150) GeV (M y M 1 )=(400, 200) GeV 




FIG. 8: Total cross section for the amplitude (+ -) as a function of the dimensionless parameter 
5ll and ^/s 77 = 410 GeV. The values of the others parameters are given in the legends. 



6.8 x 10~ 8 3.9 x 3.9 x 10~ 7 while the black triangle-up shaped points (red in color), that 
satisfy T(/i — > ej) < 1.2 x 1CT 11 cover a more restricted part. The grey circle shaped 
points (magenta in color) are determined imposing the condition that the total cross section 
multiplied by the luminosity gives more than five events per year. We can note two things: 
the signal's points overlap with the "/•«, e" region only on the tail of the red region extending 



to higher values o 
discussed in Ref. 



5ll- This tail is due to some peculiar cancellation between diagrams, as 
9^, thus we can say that this possible final state is almost excluded. The 
/i, t and the e, r final state are not excluded but they generally require a high-mass splitting 
Sll > 10- 1 . 

From the point of view of the supersymmetric seesaw mechanism described in Section [HJ 
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FIG. 9: Scatter plot in the plane (5ll, rfi) of: (a) the experimental bounds from [i — > and 
r — > ^7 (allowed regions with circular dots); (b) regions where the signal can give at least five 
events at year for two sets of gaugino masses. The energy is ^/s 77 = 128 GeV and the luminosity 
L = 136 fb^ 1 yr" 1 . 



these values can be realized in nature only under some restricted conditions [2J| : the matrix 
Y v from the seesaw mechanism neutrino masses and mixing is ambiguous up to a complex, 
orthogonal matrix R ^j]. Usually this matrix is taken to be real or identical to the unit 
matrix. In the case of a quasi-degenerate neutrino mass spectrum, R being complex allows 
for values of Am 2 larger by 5-8 orders of magnitude relative to the case of R being real or 
the unit matrix |7j. 
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2E =410 GeV 
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FIG. 10: Scatter plot in the plane (8LL,fh) of: (a) the experimental bounds from fi —* ej and 
r — > /i7 (allowed regions with circular dots); (b) regions where the signal can give at least five 
events at year for two sets of gaugino masses. The energy is ^/s 77 = 410 GeV and the luminosity 
L = 341 ft)" 1 yr _1 . 

V. STANDARD MODEL BACKGROUND 

Production of charged leptons will be copious in 77 collisions, and the SM provides several 
processes that can mimic er, fir final states. Let us see how to reduce the most important 
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TABLE I: Total cross section without and with cuts described in the text for the background 



processes Eqs. (|24l26j) . 



2E (GeV) 


77 — > tt — ► revu 


77 — > WW — * ervu 


77 — ► Tree 


200 


a 


0.58 fb 


2.3xl0 _1 


36.7 pb 




@ cut 


1.49xl0~ 6 fb 


// 


4.4xl0~ 2 fb 


300 


a 


3.1 fb 


0.48 pb 


38.9 pb 






16.3xl0~ 6 fb 


// 


3.7xl0~ 2 fb 


400 


a 


4.9 fb 


0.69 pb 


39.5 pb 






3.9xl0~ 4 fb 


2.1 x 10~ 2 fb 


2.9xl0~ 2 fb 


500 


a 


6.1 fb 


0.77 pb 


39.9 pb 




&cut 


9.7xl0~ 4 fb 


1 x 10" 1 fb 


2.4xl0~ 2 fb 



contributions that are : 

(a) 77 — ► t~t + —> T~v e i? T e + (24) 

(6) 77 -> W-*W + * -> t~v t e + u e (25) 

(c) 77 — > e + e~T + r~ (26) 

with similar processes for the production of \xr pairs. As we have seen, the e/x final state, 
which is the easiest to reconstruct from the experimental point of view, is almost completely 
excluded by the strong bounds from the non observation of the radiative decay /i — > ej. 
Thus we are bound to consider signals with a tau in the final state. Taus can, in principle, 
be reconstructed looking at the associated leptonic decay r — > Ivv and at the hadronic 
decay r — > h^-r® . The cross sections of processes in Eqs. (I24H25)) depend on initial photon 
polarizations, while the reaction in Eq. (J2*rjj) is almost insensitive to photon helicities. We 
use the program COMPHEP |2^, and the CB spectra with z > 0.8z max . In Table U] we 
give the values of the cross sections after the application of kinematical cuts (contributions 
of the charge-conjugate processes are also included). Tau pair production, W gauge bosons 
pair production and four charged fermion production are known to have very large cross 
sections, at the level of hundreds of picobarn at the CM energy of ILC, orders of magnitude 
larger than the signal in the most favourable regions of the parameter space. However the 
signal is characterized by two back-to-back leptons with the energy of the beams without 
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missing transverse momentum and energy. These characteristics provide also indications on 
the necessary kinematical cuts to be applied to the background processes. 

The helicity amplitudes which dominate the signal (+, — ) and (— , +), are peaked along 
the collision axis. Most of the background is also concentrated in this region. So we apply 
the angular cut | cos(0)| < 0.9 (6 < 25.8°) both to the signal and to the background. We also 
impose the back-to-back condition on the background processes, requiring 180° — 9w < 5°. 
Using in addition the condition that one of the event hemisphere should consist of a single 
muon or electron with energy close to E 7 , final leptons are required to have energy at least 
85% of the maximum photon energy E^ ax = y max E . 

As can be seen from Table HI after these cuts are applied, process (a) is suppressed 
because tau pairs are almost produced along the collision axis, and process (b) is completely 
excluded, at least for low energies, because the leptons from the decay of W are less energetic 
and cannot survive to the energy cut. Instead, due to the well known rapid growth of the 
77 — > WW cross section above threshold, at 400 and 500 GeV CMF energy, these cuts are 
not enough to suppress the background, giving cross sections of 2.1 x 10~ 2 fb and 10 _1 fb 
respectively. However with a cut on the transverse momentum of the electron > 15 GeV 
the cross section, at a/s^ = 500 GeV, is reduced to 2.1 x 10~ 2 fb, while for pf, > 20 GeV, 
the contribution is eliminated. 

Reaction (c) turns out to provide the most dangerous background. In this case the 
results were obtained with a MonteCarlo code developed by some of the authors [^J, which 
uses some compact analytical expressions for the diagrams with the exchange of space-like 
photons. The configuration that mimics the signal arises if one er pair is emitted at small 
angle with respect the collision axis and is not detected (we require 

^untagged < 

the other pair is tagged. This configuration is determined by multi-peripheral diagrams as 
a consequence of t-channel poles at small angles. The detected pair presents characteristics 
very similar to those of the signal, and though the cross section is effectively reduced by 
orders of magnitudes, it is still at the level of 10~ 2 fb, thus remaining competitive with the 
signal cross section. However, at a final step, this background can be estimated from the 
data by requiring instead that the detected tau and electron be of the same charge, and 
eventually subtracted. After the cuts discussed above this is the only significant background 
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contribution which remains. We consider the statistical significance 

La" 



SS 



cut 



v 7 ^ 



BG 
cut 



and requiring SS > 3 we obtain af ut > 5.4 x KT 2 fb at v /s^ = 200 GeV and a s cut > 2.5 x lO" 2 
fb at ^J~Se~e — 500 GeV, using the simulated annual luminosity for TESLA. By inspection of 
Fig. [Hand Fig.|E]it is seen that this condition is in both cases satisfied if 5ll > 10 _1 with the 
values of the other SUSY parameters as specified previously. This region of the parameter 
space is allowed for the er, fir channels as can be seen in Fig. El and Fig. E3 



VI. SUMMARY AND CONCLUSIONS 



We have studied the lepton flavor violating reactions 77 — > ££' (£, £' — e, //, r, £ 7^ £') 
which arise at the one loop order of perturbation theory and which will be of interest for the 
77 option of the future ILC The LFV mechanism is provided by low energy R-conserving 
supersymmetry with non diagonal slepton mass matrices. The origin of the non diagonal 
entries of the charged slepton mass matrices can be ascribed to a SUSY seesaw mechanism 
with mSugra boundary conditions, a theoretical scenario that has attracted much attention 
in the literature in recent years. We have studied the signal in a model independent way in 
order to pin down regions of the SUSY parameter space, the (m^, Sll) plane, allowed by the 
present experimental limits. 

We have shown that in the range 200-500 GeV for the center of mass energy of the basic 
electron collider that produces photon beams, the cross section of the signal is (7(77 — > 
££') = (9(10 _1 — 10~ 2 ) fb for sparticle masses in the range 90 — 200 GeV that correspond to 
a light SUSY spectrum somehow hinted to by fits on standard model parameters and SUSY 
benchmark points. Observation at a PC of 77 — > er, (fir) is not excluded by present bounds 
on the radiative lepton decays r — > e7, r — > /i7 which do not constrain the parameter space 
strongly enough, nonetheless a relative mass splitting Sll = Am 2 /m 2 at least of order 10" 1 is 
required, a value that can be obtained in the SUSY seesaw framework but only within some 
particular model. The efi final state is almost excluded because of the stronger constraint 
provided by the upper bound on the branching ratio given by the non observation of fi — » ey 
which is four orders of magnitude smaller than those provided by the non observation of 
t — > ej, t — > ^7. 
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The nice signal's feature of having two back-to-back high energy leptons in the final 
state can be somewhat altered by the non-mono chromaticity of the photon beams produced 
via Compton backscattering. So we have restricted the numerical analysis to the high- 
energy part of the luminosity spectrum of the photon collider: on the other hand this part 
corresponds to collisions of almost monochromatic and polarized photons, with an integrated 
luminosity close to that of the basic lepton collider. Under the very same conditions we have 
studied the standard model background and shown that with suitable cuts it can be taken 
at the level of aback ~ O(10~ 2 ) fb. The process 77 — > ee{\i\i)rr presents a configuration with 
an undetected er pair emitted at small angle along the collision axis and with the detected 
pair of high energy leptons almost back-to-back, has a potentially large cross section which 
can easily mimic the LFV signal. We have considered the signal statistical significance and 
found that one can obtain SS > 3 provided that 5ll ^ 10 _1 . 
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APPENDIX A: HELICITY AMPLITUDES 

In this appendix we present explicit expressions for the helicity amplitudes of the diagrams 
depicted in Fig. [TJ 

(a) Penguin diagrams 

These are the diagrams depicted in part (a) of Fig. ^ We have two types of contribu- 
tions: the chargino-sneutrino loop and the slepton-neutralino loop: 

1. Chargino-sneutrino 

Mf = +i e -^-X* x ' (. = 1,2) (Al) 
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xr 

xr 



+ sinr [(C M| - 2C 00 ) + s(d + C n + C 12 )] 



+ sinr (C M|-2C 00 ) 
-sin0* (C M| - 2C 00 ) 
+ sinr (C M|-2C 00 ) 



i - cos e* 

1 + cose* 



xt 



X,— = 



= +sinr(C M|-2Coo) 



+ sinr (C M|-2C o) 
-sin0* (C M|-2C 00 ) 



1 — COS ( 



1 + cos 9* 

+ sine* [(C M| - 2C 00 ) + s(C 2 + C 12 + C 22 )] 



In the above expressions the three-point form factors (C's) are to be evaluated 
with the following arguments: 

C... = CJt,0,0,ml o ,m 2 7 ,m 2 7 ) (i = 1) 



C... = C..(0, 0, t, m£ , mj, mf) (i = 2) 



2. Slepton-neutralino 



(A2) 



xr 
xr 

xr 

xr 
xr 
xr 

xr 



+ sine 

+ sine* |+Coo 
+ sine* 
+ sine* 



+C 00 + " cos^Xd + C n + C12) 



1 - cose* 



(l-cos0*)(Ci + Cn + C7i2) 



1 + cos e* 

-Coo - |(1 + cose*)(C 1 + C n + C 12 )] 
+C 00 - |(1 + cose*)(C 1 + C n + C 12 )] 



+ sine 

+ sine* |+Coo 

+ sine 
+ sine 



+C 00 - -(1 + cose*)(C 2 + C 12 + C 22 ) 



1 - cose* 



-(i-co S e*)(c 2 + c 12 + c 22 ) 



1 + cos e 

-Coo + |(1 + cose*)(C 2 + C 12 + C 22 )] 

+c 00 + |(i - co S e*)(c 2 + Ci 2 + C22)] 
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In the above expressions the three-point form factors (C's) are to be evaluated 
with the following arguments: 



a = r. 



C... = C...(0,0,t,m|,M|,M|) (i = 2) 



(b) Self-energy diagrams 

1. external leg corrections 

• slepton-neutralino 

(M 3a + M ia ) xx ' = 

• chargino-sneutrino 

(M 3b + M 4b ) xy 



167T 



O^Ol q 2 [B + B 1 ] (0, mf, mfo) X xx ' (A3) 



16vr 5 



(0^2 [B + SJ (0,mlm 2 ~)X xx ' 



2. t-channel correction .M5 = M.§ a + A^sj, 
• slepton-neutralino 



M 



5 a 



16tt 5 



The helicity factor is the same in this case: 



X 



sin 8* 
1 + cos 9* 



1 + AA' + A - A' 



AA' cos6* 



x ++ 


= +sin#* 


x + - 


= -hsin^* 


x- + 


= -sin0* 


x— 


= +sin6>* 



1 - cos e* 

1 + cos 9* 



(c) Sea- Gull and Box diagrams 

These are the diagrams depicted in part (c) of Fig. [TJ 
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(A4) 



O e ~oOio \B (t, mf, m|o) + B 1 (t, mf, m| )l X xx ' (A5) 



chargino-sneutrino 

M xx ' = -i^{0ff [Bo{t,mlm%) + B^mlm^)} X xx ' (A6) 



(Seagull) The seagull type diagram has only a contribution from a slepton- 
neutralino loop (.Me). 



Mf = 



(Scalar Box) 

This is the slepton neutralino box diagram in part(c) of Fig. Q {M.7) 



M 



xx 
7 



16n'- 



l + AA' 



002 



A - A' 



AA' cos 6* ( A)o + Arai + ^002 + D, 



003/ 



— o AA sin" 6* (D112 + 2,0x22 + 2D 123 + £222 + 2D 223 + D 



233 



+D 2 + 2{D 12 + D 22 D 



23 J 



(A7) 



The D form factors appearing in the above formula for the scalar box diagram 
are to be evaluated with the following arguments: 

D... = D...(0,0,0,0,i, s,m 2 x o,m~,m~,m~) 
Chargino-sneutrino loop 

This is the box diagram involving fermions depicted as (Ads) i n part(c) of Fig. 



M 



XX' 



W\2 



. i ij ijk 



(A8) 



The various (non-zero) coefficients multilplying the four-point loop form factors 
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(Di, Dij, Dijk) are given below: 

(0) AV = -M| | sin 0* (1 + AA' + A - A' + 2AA' cos 0*) 

(1) AA ' = +M^A'sin0*(l - Acos0*) 

(2) AA ' = +M| | sin 0* (1 + AA' - (A - A') - 2AA' cos 0*) 

(3) AV = -M||Asin0* (1 + A'cos0*) 
(00) AA ' = +^ sin0* [1 + AA' + (A - A') + 2AA'cos0*] 

(12) AA ' = +^ sin0*A'(l- A)(l + cos0*) 
s 2 

(13) AA ' = -— sin0* [1 + AA' + (A - A') + 2AA'cos0*] 

(22) AA ' = --(A- A')(l- A) sin0*(l + cos0*) 
8 

s 2 

(23) AA ' = -— A(l + A') sin0*(l + cos0*) 

(001) AA ' = -ssin0*[A + A' + A'(l- Acos0*)] 

(002) AA ' = (|sin0*)[-(l + AA') + A(l + A'cos0*) - A'(l - Acos0*)] 

(003) AA ' = +s sin 0* [A(l + A' cos 0*) + A + A'] 

s 2 

(112) AA ' = + — A'(l-A)sin0*(l + cos0*) 
s 2 

(113) AA ' = +_A'sin0*(l - Acos0*) 
2 

s 2 

(133) AA ' = -— Asin0*(l + A'cos0*) 

Z 

(122) AA ' = -— (3A'- A)(l- A)sin0*(l + cos0*) 

A'(l - A cos0*) - \(\ + AA' + A + A') 

(132) AA ' = 5 X+ '^ sin0*(l + cos0*) 
z 

s 2 

(222) AA ' = —-- (A - A')(l - A) sin0*(l + cos0*) 

8 

s 2 

(223) AA ' = - — A(l + A')(3- A) sin0*(l + cos 0*) 

8 

s 2 

(233) AA ' = -— A(l + A') sin0*(l + cos0*) 
The D form factors appearing in the above formula for the chargino sneutrino box 



(123) 



AA' 



+ — sin0* 
z 



25 



diagram are to be evaluated with the following arguments: 

D .. = A..(0, 0, 0, 0, t, s, m|, m|, m%, m 2 ~) 

A diagram with a LFV and a LFC scalar line, for example, is described by the propagators 
of Eq. (jSJ) and Eq. so that the loop coefficients in the amplitudes are a sum of four 
integrals, while in the diagrams with only a LFV line, are a sum of two. The scalar two 
point function Bo and the tensor coefficients B\, Coo that appear in the electroweak penguins 
are ultra-violet divergent, but the amplitudes are finite due the ortogonality of the slepton 
mixing matrix. The rule for obtaining the helicity amplitudes for the exchanged diagrams 
from those of the direct diagrams: 

MZ h .(^ e\ cos e*) = M x 2 ect (-^o*,- cose*) (as) 

The loop form factors are exchanged accordingly to the same rule: (cos 6* — > — cos 8*, and 
sin6>* — > — sin6>*). 
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